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Abstract: We consider the hydrodynamics of supersymmetric fluids. Supersymmetry is 
broken spontaneously and the low energy spectrum includes a fermionic massless mode, 
the phonino. We use two complementary approaches to describe the system: First, we 
construct a generating functional from which we derive the equations of motion of the 
fluid and of the phonino propagating through the fluid. We write the form of the leading 
corrections in the derivative expansion, and show that the so called diffusion terms in the 
supercurrent are in fact not dissipative. Second, we use an effective field theory approach 
which utilizes a non-linear realization of supersymmetry to analyze the interactions between 
phoninos and phonons, and demonstrate the conservation of entropy in ideal fluids. We 
comment on possible phenomenological consequences for gravitino physics in the early 
universe. 
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1. Introduction 

Spontaneously broken symmetries play an important role in the physics of fluids. Phonons 
can be considered as the Goldstone modes of Lorentz symmetry, which is broken by the 
choice of the fluid's rest frame. Another example is the phenomenon of superfluidity, 
where the spontaneous breaking of a £7(1) global symmetry leads to the appearance of a 
non-dissipative flow in the fluid. 

Supersymmetry is broken at finite temperature since the thermal ensemble differen- 
tiates between the statistics of fermions and bosons. Interestingly, this symmetry is not 
restored at high temperature |lj. It is well known that when supersymmetry is sponta- 
neously broken, the low energy spectrum contains a fermionic massless mode, known at 
zero temperature as the goldstino. A similar situation happens in super symmetric fluids at 
finite temperature, where a long range fermionic fluctuation appears. This mode is called 
phonino. The existence of this mode was demonstrated by using the real-time formalism 
in |2j. Its dispersion relation was calculated for strongly coupled supersymmetric plasma 
using the AdS/CFT correspondence in ||]. At leading order in momenta, the phonino 
satisfies a linear dispersion relation with velocity (see also pi M): 

(1.1) 
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where P is the pressure and e is the energy density. In simple cases where P oc e, the speed 
of the phonino is the square of the speed of (first) sound. At zero temperature |P| = |e| 
and the phonino is the goldstino, a free relativistic fermion. 

The existence of this mode can be understood as a consequence of a super symmetric 
Ward-Takahashi identity for the supercurrent two-point function || : 

d, (T{S»(x)S»{y)}) = 5^(x - y) (T up ) a p . (1.2) 

Going to momentum space and assuming a constant energy-momentum tensor the identity 
becomes 

-ik^ s = (T»») a p . (1.3) 

In order to satisfy this identity for all momenta we need that T$s ~ which is possible if 
there is a massless fermionic mode, the goldstino/phonino. Let us assume that the system 
is at thermal equilibrium and that the energy-momentum tensor takes the form of that of 
an ideal fluid 

(Tn = dmg(s,P,P,P) . (1.4) 

As in the case of spontaneous breaking of supersymmetry, we expand the supercurrent 
(S p ) at low energies in derivatives of the phonino field G (G) : 

S° = iF t a°G + --- , S i = iF s a i G + --- , (1.5) 

where we have assumed rotational invariance but not Lorentz invariance. In these conven- 



tions G is a Weyl fermion. The identity (1.3) becomes 



- i(F t k a° + Fsha 1 ) (GG) F t a° = ea°, (1.6) 

- i(F t k a° + F.kia*) (GG) F s a l = -Pa 1 , (1.7) 

The equations require that the propagator of the phonino has the form 

{GG) ~ F t kl - vlk} ' (L8) 

where 

F t = e, F s = -P, v G = -£ = —. (1.9) 

There is some arbitrariness in the normalization, this choice corresponds to the definition 
of the supercurrent 

S» = iT»a v G + --- (1.10) 
such that the conservation of the supercurrent becomes 

= 3^ = iI>^G + --- (1.11) 



1 In our conventions the supercurrent has a relative factor l/y/2 respect to the one originally used. This 
removes a factor of 2 from the right hand side of the WT identity. 
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At zero temperature Ty = — e5y and the propagator of the phonino becomes that of 
the goldstino. In this case the energy density is determined by the scale of supersymmetry 
breaking e = |/| 2 . When Lorentz invariance is broken, \vq\ 7^ 1 is the velocity of the 
phonino. Note also that the conformal Ward identity for the supercurrent agrees with the 
tracelessness of the energy-momentum tensor in a CFT 

a^oc(T^G+--- = . (1.12) 

These two equations are valid to lowest order at the ideal level even out of equilibrium, 
although there can be derivative corrections beyond the ideal level. 

In this work we will be studying the hydrodynamic regime of supersymmetric field 
theories, that is the the combined effective theory of the phonino and the fluid. In particular 
we will be interested in the motion of the phonino through the fluid itself and in its 
interactions with other long range modes in the fluid, the phonons. In order to describe 
the motion of the phonino through the fluid, we will begin by constructing the effective 
action following a method inspired by the recent works || [l(| 2 to derive relations between 
transport coefficients from a generating functional depending on external metric and gauge 
fields. We will adapt the method to include gravitino sources and the phonino. From the 
resulting generating functional one can derive the equations of motion of the fluid and of 
the phonino propagating through the fluid. Since the gravitino sources are included, it can 
also be directly interpreted as the contribution to the effective action of the gravitino from 
the fluid. As we will see, a mass for the gravitino appears at the order of the first derivative 
corrections of the fluid. 



We note that supersymmetric hydrodynamics has been studied in [13]. Our analysis 



is rather different. In [ 1 3 1 the concepts of "classical" supersymmetric charge and chemical 
potential are introduced and are used to define the constitutive relations with which one 
can express the supercurrent. This is done is analogy to the hydrodynamic description 
of bosonic currents such as the stress-energy tensor and global symmetry currents. How- 
ever, one cannot define consistently classical fermionic charges as the expectation values of 
fermionic operators is zero. We therefore propose that the correct analysis of the hydro- 
dynamics of supersymmetric field theories is not that of a flow of a fermionic charge, but 
rather the dynamics of the phonino in the normal fluid. 

Although in principle one can derive all the relevant information from this generating 
functional, in order to study interactions of the phonino with the phonons in the fluid 
it is better to use a slightly different set of variables, that in particular will allow us to 
make contact with the effective actions at zero temperature of the goldstino existing in the 
literature. The approach is based on the Akulov-Volkov formalism for non-linearly realized 
supersymmetry (NLRS) [0|. The main results of this analysis are the following: first, we 
introduce a formulation to relate NLRS to hydrodynamics via conserved currents which is 
shown to be consistent on both ends - it is a valid hydrodynamics description of ideal fluids 
in the sense that entropy is shown to be conserved, and it has the correct transformation 
properties under supersymmetry. The important new constraint we impose is that both 



2 See also pi], where the partition function is constructed from a microscopic theory with free fermions. 
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the currents and their conservation equations transform in a way which is consistent with 
the SUSY algebra. This leads to a unique construction of the physical currents in terms of 
currents dressed with appropriate goldstino factors. Since this construction is dictated by 
symmetry, it is valid even in the absence of a Lagrangian description for the microscopic 
system. Secondly, we use this formalism to study the dynamics of phoninos in ideal fluids. 
We find the phonino dispersion relations and compute the expressions for the leading 
interactions between phoninos and phonons. 

The paper is divided into two sections - one for each of the two formalisms mentioned 
above. We conclude with a short discussion of the phenomenological implications of this 
analysis. A review of non-linearly realized supersymmetry and details of some of the 
computations appearing in the paper are given in the appendix. 



2. Generating functional with a phonino 

We will consider a relativistic fluid in the presence of external sources. We will assume 
that all the microscopic degrees of freedom have been integrated out and that we can use 
a hydrodynamic approximation to describe the properties of the system. If the underly- 
ing theory is super symmetric we expect that a collective massless fermionic fluctuation 
exists, the phonino. The motion of the phonino should be included in the hydrodynamic 
description, and as a first approximation we can treat the phonino as propagating through 
a thermal medium but not affecting its properties. The effective action of the phonino 
will depend on the properties of the medium and on external sources, and we can use 
supersymmetry to constrain its form. The partition function of the theory will take the 
form 

Z= I DGDGe iW[G] , (2.1) 



where G is the phonino and W is determined by the temperature and possible external 
sources. In principle W is a non-local functional, but it has a local expansion when static 
configurations (in the rest frame) are considered as explained in |9[ 1C]. 



In order to describe the effect of the medium, the temperature and the chemical po- 
tential are treated as external sources and W becomes the generating functional for them. 
In the absence of chemical potential, the generating functional takes the form 

W = f d d x^g~P{T) + ■■■ (2.2) 



Where T = Tq/^/—V 2 = Tq/ ^—g^V^V^ is the temperature and V M is proportional to 
the velocity of the fluid. 

The energy-momentum tensor is 

2 5W dP V» V v 

= = P<T + T— —=—= = PgT + TsuW. (2.3) 

V=5 og^u dT y/-V 2 V-V 2 

The equations of motion can be deduced from diffeomorphism invariance of the underlying 
theory. Although the background metric breaks physical diffeomorphism invariance, the 
partition function is invariant up to anomalies under a combined transformation of the 
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dynamical fields and the external sources. This implies that the generating functional 
should be invariant up to anomalies under transformations of the sources 



S a 9fiu = y^a v + V^. (2.4) 

Then, 

= 8 a W = / d d x- S a g„ v => V M r^ = 0. (2.5) 

Suppose now that there is a global symmetry with a chemical potential 3 

M = ^- (2.6) 

The current is defined as 

1 5W dP u . . 

JM = ; f a = ~ ; = /°^- ( 2 - 7 ) 

invariance of the generating functional under gauge transformations bxA^ = d^X imposes 
the following condition 

= 5 X W = - j dtxy/^gXVpJ 1 *. (2.8) 

If the global symmetry is spontaneously broken there is a Goldstone boson eft. Under a 
global symmetry transformation 6\(f) = A. The chemical potential can be defined as 

,- V(V^« ( , 9) 

which is explicitly gauge invariant and gives the Josephson condition u^V ^(ft = —\x in the 
absence of sources. A new scalar quantity on which the generating functional can depend 
is (Ap - V M 0) 2 . 

If the underlying theory is supersymmetric, we can extend the same arguments to 
supersymmetric transformations. The difference is that now we should consider more 
general external sources, including a gravitino. The generating functional is a function of 
the vierbein, spin connection and gravitino background values (we do not consider global 
bosonic symmetries for the moment) 

W[e%u?,*«,*% = J d d xV^P(T,fi a ,fin, (2.10) 



Here a refers to spinor indices and we denote 

The generating functional should be understood as having an expansion in the Grassmann 
fields. ^ to be a Weyl spinor in the notation of 14]. 



3 This expression is not explicitly gauge- invariant. A gauge-invariant expression involves a Wilson line 
along the time direction defined by the Killing vector V^d^. Fixing the gauge so that only time- independent 
transformations are allowed, the chemical potential takes the form (|2.q). 
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Under a super symmetric transformation, the variation of the sources is (suppressing 
spinor indices) 



= 5^W -- 
Therefore, the conservation 



<%e£ = z(vIV^-^<%), 

^ = -22^, 



(2.12) 
(2.13) 
(2.14) 
(2.15) 



The derivative includes the spin connection. These transformations are actually special 
for Poincare symmetry. Had we considered for instance a theory in anti-de Sitter space the 
transformations would be different. Introducing external gauge fields will also modify the 
supersymmetric transformations. 

A general supersymmetric transformations of the generating functional leads to 



-I 



d 4 x 



SW . a 5W .= rT SW 



(2.16) 



= J (fx [iT^^ + 2^5" + h.c] 
equations for the supercurrent are: 



(2.17) 



2V fX S» = iT^o^v, 2V fl S fM = -iT^V^ . (2.18) 

Now we can take into account the effective spontaneous breaking of supersymmetry, as 
it happens at finite temperature, by introducing the phonino field, whose transformation 
under supersymmetry is 

6zG = S, 5 ( G = l (2.19) 

The scale of supersymmetry breaking appears in the definition of the supercurrent. This 
field is analogous to the phase of the condensate in spontaneous breaking of global sym- 
metries. 

The effective action of the phonino will be determined by the original generating func- 
tional in an analogous way as some terms in the effective action of Goldstone bosons are 
determined by the generating functional for gauge fields. With the phonino we can con- 
struct the SUSY-invariant combinations 

^ = + 2X>„G, $ lt = % + ZD ll G, (2.20) 

for the gravitino sources. For the vierbein we can form the combination 

E£ = e%- i(*^a a G - G<r<%) - I^V^G - Ga a V^G). (2.21) 

The transformation of the dressed vierbein is 

8<:E? = -2i{V^Ga a i - fr a V^G). (2.22) 
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We can also write it as 

E^ = e^-i(^a a G-Ga a ^) . (2.23) 

This choice for the dressing of the vierbein leads to an action whose supersymmetry vari- 
ation is independent of the gravitino and, as shown below, gives the correct equations of 
motion. This approach is similar to the AV formalism we will use in the next sections, but 
the hydrodynamic variables will be described differently. Note also, that in the current 
approach the SUSY transformations are linear and the variation of the generating func- 
tional gives the conservations laws, while in the AV formalism the SUSY transformations 
are nonlinear under which the effective action is invariant. 

The dependence on the metric is changed to a dependence on the covariant combination 
9\w — > G^ v = ri^E^E^. Then, the pressure depends on the dressed temperature 

(2.24) 



where 

G P W = g^VV (l - 2iu a (fia a G - Ga a j2) - Vab (fia a G - Ga a j2)(fia b G - Ga b p)) 

(2.25) 

The velocity and the fermionic chemical potential are defined in the usual way, as we will 
also do with the temperature 

^ = -=, n = vPij> lil T=^=. (2.26) 
yf-V 2 yf-V 2 

where we have defined V 2 = guuV^V^. 

In the absence of sources the definition of the chemical potential gives a relation u^ipn = 
u^V^G = fx similar to the one found in a superfluid. In the generating functional the 
covariant form of the chemical potential should appear 

The pressure is a function of the temperature and fermionic bilinears P(T,i() 2 , ft, fi), 
but not all the terms are of the same order. We are doing an expansion in small derivatives 
V M ~ e. In order to keep all contributions to the vierbein of the same order ~ 1 and 
assuming that the velocities are of the same order ~ 1, we should fix the scaling of the 
phonino to G ~ 1 /y/c. In this case Vv ~ so bilinear contributions to the stress-energy 
tensor are suppressed ~ e. Similarly, this scaling behavior ensures that contributions to the 
supercurrent which are linear in tp are suppressed by a factor of e with respect to the leading 
contribution, which is proportional to G. We conclude that the contributions related to 
tp and fi are of order e, the same as the contributions of bosonic fields that contain one 
derivative. 

We now define the generating functional to leading order as 

W = j d 4 xEP(f), (2.28) 
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where E is the veilbein determinant. Under a supersymmetric variation 

5^W = J d A xE {2iT^ a V^Ga a i + h.c) , (2.29) 

where the energy-momentum tensor has been defined as 

1 SW 

Tn = (2.30) 

a E5E a a y ' 

Therefore, supersymmetry requires that the following condition is satisfied 

T^ a V^Ga a = . (2.31) 

Using the condition ( 2.31| ) and the Fierz identities 

X° a $ = -^cr a x, (2.32) 

the generating functional satisfies the Wess-Zumino consistency condition 

(6 n 5 e - 5^ V )W = 2i J d A x ET^ a V^{r,a a i - £a a f}). (2.33) 

Since transforms under diffeomorphisms in the same way as e^, diffeomorphism 
invariance of the action requires that T£ is conserved. The conservation equation is, in flat 
space 

d ll (ET» a ) = 0. (2.34) 

Although there is a factor depending on the vierbein, its derivative is zero in the absence 
of gravitino sources 

EE^T^El = -AiEE» b T» a d^Ga a d v G = . (2.35) 
Where we have used that 

T il a d VL Go a = . (2.36) 

We can use the vierbeins to define the energy-momentum purely in the orthogonal 
frame 

T% = E%T» h . (2.37) 
The explicit form of the energy-momentum tensor is 

We can simplify further this expression using that 



E ayn = ya_ iV H^ _ GV^) = V° + yf-V^C = \f-V*{u a + C°), (2.39) 



where V a = efY^ as usual and we have defined 



C* = -i(fia a G - Ga a Jl) . (2.40) 
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Note that the dressed temperature is 

f = T . (2.41) 

Assuming there is no external gravitino, there are only four independent components 
in the Grassmann field G, so there can be up to 0(C 2 ) terms. For the temperature this 
means 

f = T (l + u a ( a + i W + 3uV)CaO>) . (2.42) 

The pressure then can also be expanded as 

r) P ( 1 \ 1 f) 2 P 

P(T) = P(T) + — \v?Ca + z(v ab + 3u°« 6 )CaJ + -^u a u b U b . (2.43) 

Finally, the energy momentum tensor takes the form 

T ab = Vab (P(T) + amXc + (a 2 v cd + a 3 u c u d )( c ( d ) 
+ u a u b (t^ + hu c ( c + (b 2 r] cd + hu c u d )CcU 

+ Mb + U b (a) + Cl U%)j + CaCbT^. (2.44) 

Where 

dP a\ dai 

ai= 9T' a2 = Y' «3 = 3a 1 + — , 

bi=T { 2ai+ w)> b2= r 63=T r i+ 29T + ^J' 

ci = 6i. (2.45) 

The first derivative corrections depending on the phonino are easy to compute. First 
note that the field tp^ is in a (|, 0) <8> (|, |) = (0, ^) © (l, ^) representation of the Lorentz 
group, while V> M is in (0, ^) ^) = 0) © (|, l). Therefore, there are four possible 
scalar bilinear combinations 

G >xv i>ni>v, GT$$ V , %a^ v , ^aWfa, (2.46) 

all of which are of order ~ e in our expansion. We can add these terms in the pressure 
with arbitrary coefficients depending on the dressed temperature 

P = P {f) + Pi(f )A 2 + P 2 (f)G^^„ + P 3 (f)^a^v + h.c. + ■■■ (2.47) 

Note, that the new terms can also be interpreted as bilinear terms for the gravitino field 
induced by the fluid. In particular, the term proportional to Ps(T) has the form of a 
mass for a Rarita-Schwinger field, while P 2 (T) and could be interpreted as a mass for the 
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spin-1/2 components of the gravitino field. Note that both coefficients can be complex in 
general, for instance 

P 2 {f)G^^ v +h.c. = ReP 2 (f) {CT%i> v + G^tP^)+iImP 2 {f) [G^^^ v - GT^u) ■ 

(2.48) 

Both terms can contribute to the dispersion relation of the phonino, adding a contribution 
quadratic in the momentum which has a real part proportional to Re P 2 (T) and a imaginary 
part proportional to Im P 2 (T). In the conformal case P 2 (T) and i"3(T) are not independent. 
P 2 (T) has been calculated for particular examples in the AdS / CFT framework ||] and 
it modifies the phonino dispersion relation. It was found that for N = 4 SYM, Re P 2 (T) = 
and Im P 2 (T) 7^ 0, implying that this is a damped mode. However, it is clear from our 
construction that although it apppears at the first derivative order it is not a dissipative 
term in the sense that it does not contribute to the entropy current. 

In the following sections we will introduce a slightly different formalism that we will 
use to compute the interactions between phoninos and phonons. Although the results 
are equivalent in the two formalisms and there are many parallelisms, the differences are 
important enough that we will give a detailed description of it. 



3. Effective field theory approach for phonino dynamics 

Phonino dynamics can be studied using an effective field theory. This approach is based 
on a Lagrangian description, and therefore valid only for ideal fluids, but it will give a 
new perspective on the results of the generating functional formalism and enable us to 
investigate the interactions with the other long range modes in the fluid, the phonons. 
As in the previous section, we write the most general action allowed by symmetry, and 
the effects of spontaneously broken supersymmetry are manifested in the appearance of 
the fermionic zero mode. The framework we will use to incorporate the phonino into 
the dynamics of the fluid is the Akulov-Volkov (AV) formalism for non-linearly realized 
supersymmetry (NLRS), which is reviewed in appendix [A| |l2|]. 4 At this point we would 
like to emphasize three points related to this framework. First, we define the standard 
non-linear realization of supersymmetry for a generic field <p as follows: 

Stf = -i {Ga^l - ia^G) drf . (3.1) 

In the appendix it is shown that the transformation properties of the Goldstino are such 
that this transformation rule satisfies the supersymmetry algebra. 
Secondly, we introduce the AV vierbein 

\ a = ^ + i {Ga^G - d^Ga a G) . (3.2) 



This expression is essentially equivalent to the construction discussed in (2.23) but we will 



use this notation in this section in order to avoid confusion with the standard literature in 



4 Another possible formalism for NLRS is the framework of constrained superfields which has recently 
gained renewed attention thanks to the work of E3|. We chose to avoid the usage of superfields because at 
nonzero temperature they require non-trivial boundary conditions in superspace. 
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the field. 5 This vierbein is used to define the AV covariant derivative 

V a = (A- 1 )^ , (3.3) 

which transforms according to the standard NLRS. Another well known result is that in 
order to build a manifestly SUSY invariant Lagrangian out of a non-SUSY Lagrangian, 
one has to make the following modification 

£(0,d(/)) -> A£(0,V0). (3.4) 

In other words, one has to multiply by the determinant of the AV vierbein, and replace all 
derivatives with the AV covariant derivate. 

Thirdly, we would like to discuss another aspect in which the AV formalism enters 
the low energy theory, even in the absence of a Lagrangian description of the system. As 
far as we are aware this point does not appear in the literature on the subject. Theories 
with symmetries require well defined currents, j^, which should satisfy well defined conser- 
vation equations. Since the conservation equations contain derivatives and not covariant 
derivatives, it is not trivial that the SUSY transformation of the operator d^j^ satisfies 



the SUSY algebra. In appendix [A.2| we show that in order to construct a current for which 
the SUSY transformation of both the current and its conservation equation satisfies this 
condition, one should use the following form 

f = A(A-\^f (3.5) 

where j a is an operator that in the limit where the goldstino is set to zero is equal to the 
current, and it is defined to transform according to the standard NLRS (|3.l| ). We will refer 
to such currents as "standard realization" (SR) currents, bearing in mind that these are 
merely useful notations and not necessarily physical quantities. The physical current j M 
has the same value at the zero goldstino limit, but has a different transformation law. 



In appendix A. 3 we show that this is the general form for the Noether currents in a 
theory with NLRS. However, we would like to stress that it can be used also in theories 
without a Lagrangian description, because it comes directly from constraints imposed by 
the SUSY algebra. More interestingly, in our effective field theory discussion of hydrody- 
namics we find the same form for the entropy current 

fp = A(A~ 1 )J x s a (3.6) 

and demonstrate that indeed this is a conserved current. 

3.1 A Lagrangian formulation of ideal hydrodynamics with phoninos 

In this subsection we briefly review a Lagrangian formulation of ideal hydrodynamics (see, 
e.g., |l6|, p"7|). We begin with the following Lagrangian 



C = F(B) = F ^Jdetid^d^)^ 



(3.7) 



5 Notice that in this formalism the goldstino/phonino has different transformation properties. Also, when 
using the AV formalism we will assume a flat metric. 
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where we denote B LJ = d^ 1 d^<j) J , B = ^det(B IJ ), and F(B) is some unspecified func- 
tion. F(B) is analogous to P(T) which was introduced in the generating functional for- 
malism, but it has a different physical interpretation which will be presented below. 

The (j) 1 fields can be viewed as the coordinates for the location of an element of the 
fluid in a static coordinate system. B can be written as: 

B = y/SpSf 

s » = l^euKdu^dp^d^ . (3.8) 

We note that by definition the vector s M has zero divergence. Also, it is normal to the 
gradient of the coordinate, so it has to be proportional to the velocity. Recalling that the 
velocity is normalized to u^u^ = — 1 we find the relation 

u" = ^ . (3.9) 
The canonical stress energy tensor in this system is 

T iw = _ FbB [B- 1 ) IJ d'"(t> I d u (t> J = (F - F B B)rT ~ F B Bu»u v (3.10) 
where we used the relation 

{b- 1 ) ij d^^d u <p J = A"" = rT + • (3-H) 

This relation can be understood by noting that this is a projection operator which satisfies 
A^Up = 0. We can identify the isotropic part of the tensor as the pressure and define 

P = F-F B B e = -F, (3.12) 

where e is the energy density and we denote F B = ^§ • Using the thermodynamic relations 
(in the absence of global charges) 

Ts = e + P, (3.13) 

one can identify B as the entropy density s and —F B as the temperature (up to multiplica- 
tive factors). This means that the current s M = Bu^ is in fact the entropy current, and as 
mentioned above, it is conserved by construction. This is an important point because it 
means that the system is not dissipative and therefore the description using a Lagrangian 
formalism is consistent. 

The last missing piece is the differential thermodynamic equation 

Tds = de , (3.14) 

which is necessary to show the consistency of the conservation of the entropy current when 
the equations of motions are satisfied: 

- Td^ (su 11 ) = -u^d^e ~{e + P) = u v d^ u = (3.15) 
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The final form of the action is 

S 



J d A xe(s), (3.16) 



where s = B is the entropy density and e = — F is the energy density. Note that this 
is related to the action we used in the generating functional formalism by a Legendre 
transform 

P (t - — \ - — 

\ ds J ds' 

( dP\ dP 

-i' = af) =p - T W (3 ' 17) 

This establishes the connection between the two formalisms presented in this paper. In 
both cases we define a vector which is proportional to the velocity, in the generating 
functional case and s^ 1 in the effective field theory, whose normalization is in fact the basic 
dynamical variable in the action (s = V—s 2 and T = — Sl^) and it is related via the 
Legendre transform to the other. Spontaneously broken supersymmetry will enter the 
game in both cases dressing this vector with modified vierbeins. 

The form of the Lagrangian appearing in ( |3.7| ) was dictated by the symmetries of the 
system - translations, rotations and volume conserving diffeomorphisms. Now, we would 
like to take into account NLRS as well, using the prescription given in ( |3~4 

F{B) -»• AF(B) , (3.18) 

where now 



B = ydet(V a ^V a ^) = ^f-s a s a 

§a = e abcd eijK y b(/) Iy c(f) Jy d(t) K (3 1Q) 

Using this procedure we find an action which is manifestly invariant under the non-linear 
supersymmetry transformation ( |A.1| ) and (p.l|), to be compared with the generating func- 



tional approach which is invariant only up to equations of motion under ( 2.19 ). The 
advantage of the approach taken in this section is that it contains the phonon and phonino 
as elementary fluctuations with respect to a static background, and can be used to study 
their scattering processes in this limit. 

A few comments regarding the vector s a are in place. First of all, it transforms ac- 
cording to the standard NLRS, because it is constructed using AV covariant derivatives. 
Secondly, unlike the vector defined in ([^), it is not necessarily conserved (even with 
respect to an AV covariant derivatives, because the AV covariant derivatives do not com- 
mute). Most importantly, it can be related to the vector by multiplying by AV factors: 

A(A- l )»s a = ^(A-X^A-^iA-^^A-^e^ejjKd^d^d^ = s» 

(3.20) 



- 13 - 



We will show below that the vector can be identified as the entropy current even in the 
super symmetric case. This results in two important conclusions: the entropy current in the 
super symmetric fluid is conserved by construction, and its conservation equation satisfies 
the SUSY algebra. The second point is ensured by the fact that it is written in the form 
suggested in ([?!]) 6 . 

In order to verify that the vector can indeed be identified as the entropy current, 
we derive the canonical energy-momentum tensor in (A.22), and find that it equals the SR 
stress-energy tensor dressed with AV factors 

T» a = A{A- l )ff ab = A{A- X )f (rj ab (f - F § B^j - F § Bu a u b ^j , (3.21) 
where we defined u a = £ and identified a projection operator 

IT 1 ) VV V V = A afc = rT + u b u a . (3.22) 



Making the identification as in Q3.12 ) and ( 3.13j ) for the hatted objects, we can define 
the temperature as T = —F^ and the entropy density s = B. In appendix [b] we follow 



a similar computation as in ( 3.15| ) and show that the thermodynamic relations and the 



conservation of currents indeed lead to the conservation of the entropy current also in 
the super symmetric case. 

As a side remark, we would like to mention a subtlety which is related to the definition 
of charge densities in this formalism. Unlike a current, the charge density is a frame 
dependent quantity. Using the formalism described in this section one can define two 
different normalized velocity vectors 

u ^ = — u a = 4- • (3.23) 
B B 

u a transforms according to the standard NLRS, and can be used to define a SR current 

f = pu a (3.24) 

assuming that p transforms according to the standard NLRS as well. We can also define 
the physical current using u^, but these two definitions have to be related as follows: 

f = pu» = A(A- 1 ) ( fpu a (3.25) 

Using ( 3.20Q , we find that the two definitions of charge density are related by 



p = p| . (3.26) 



6 Notice that the dressing of s a is similar to the dressing of the vector V a which was introduced in eq. 
(|2.39|) without the usage of the AV covariant derivative. 
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3.2 Phonino interactions in the fluid 



The Lagrangian formulation of classical ideal hydrodynamics can be used to find the in- 
teraction terms of the phonon with the phonino or with itself by expanding around the 
ground state: 



' + 7T 1 ) 



(3.27) 



where we identify ir 1 as the phonon, and B$ gives the value of B in the ground state. The 
leading order terms in an expansion in the ir 1 fields are: 



5C 



OF 



dB 



i.j 



o ^ 2 dBjjdB KL 



1 



A 2 



(<9/7T 



I\ 2 



5B IJ 5B KL + 

o 

+ . . . 



(3.28) 



(3.29) 



where we used ( |3.12j ) and the speed of sound is given by 

2 = OP P'(B) _ F BB B 
° s de e'(B) F B ' 

s 

Some details of this computation are given in appendix |C| The Lagrangian ( |3.28j ) can 
be used to study phonon scattering IS] and it is a good starting point to systematically 
include non-dissipative terms. 

We would like to use the same method to discuss the dynamics of the phoninos in an 
ideal fluid with NLRS, and therefore make the modification appearing in ( 3.18) ). Taking 
into account the fluctuations in the AV matrix 



SA° = i (G(7°9 M G - dfj,Ga a G) 



(3.30) 



we find new interaction terms in the expansion of the Lagrangian around the static back- 
ground: 



6 (AC) 



d(AC) 



OA, 



8A,f + 



d 2 (AC) 



dAfdBu 



SAfSBu + 



(3.31) 



The first term is the phonino kinetic term 
.8 (AC) 



Ga a d^G - 8^Ga a G) = irj b ^f ab (Ga a 8^G - 8^Ga a G) 



(3.32) 



where we used (A. 20). This kinetic term leads to the phonino dispersion relations discussed 
in the introduction. We see also that the canonically normalized phonino should be divided 

<2k. 



by a factor of y 2Too 

The next term in ( |3.31 ) gives the two-goldstinos one-phonon vertex: 



6 (AC) d iTs (\Ga°G - Ga°G) + c 2 (Ga 1 d/G - diG^G)) (djnj) 

-iTs (Ga I djG - djG^G) (3/vrj + d/vr/) . (3.33) 
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Details of the computations and computations of interaction terms with a higher number 
of phonons appear in appendix |C[ Notice that the phonino does not interact with the 
transverse mode of the phonon (or equivalently, the vortex d x ir). Using the canonically 
normalized phonino, we find that the effective coupling for the phonon- goldstino vertex 
(and similarly, the scattering of the phonino off n phonons) is proportional to Ts/e. If 
the temperature is very low in comparison with the scale of supersymmetry breaking y/J, 
the coupling will be very small as well, since e ~ |/| 2 Ts. On the other hand, at high 
temperatures the coupling will be of order one. Therefore, as the temperature is increased 
the phonino becomes more strongly coupled with the fluid. 



4. Discussion 



We presented the formalism with which the fermionic zero mode associated with the spon- 
taneous breaking of supersymmetry enters the hydrodynamics description of a fluid. The 
main ingredient was the dressing of the thermodynamics variables with a vierbein con- 
structed using the phonino/goldstino. This was demonstrated in two frameworks - the 
generating functional approach for hydrodynamics, and an effective field theory with non- 
linearly realized supersymmetry. 

The generating functional approach enabled us to write a hydrodynamics description 
for the super symmetric fluid in the presence of external sources. We considered the first 
non-dissipative derivative corrections. The supercurrent first derivative terms take the 



general form [13| 



Si = -D s ViG - D a a ij VjG . (4.1) 



The transport coefficients D s and D a correspond to are our P2 and P3 in ( 2.47| ), which 



shows that they are non-dissipative. Note that D s oc 1P2 and D a oc 1P3, so the coefficients 
can be complex in general. In a conformal theory the two transport coefficients are real 
and have the same value D s = D a . 

We have limited ourselves to the discussion of the energy-momentum tensor and the 
supercurrents, a clear extension will be to include conserved currents for global symmetries. 
In the generating functional formalism this can be implemented by including gauge fields 
and their fermionic superpartners in the supersymmetry transformations. More generally, 
the approach can be extended to larger supersymmetries, so that the sources transform as 
components of multiplets of different gauged supergravity theories. 

The effective field theory approach enabled us to write scattering vertices for the 
phoninos. Another useful result is a new prescription for dressing conserved currents with 
goldstinos in effective field theory. One simply writes the currents using objects which 
transform according to the standard realization, and multiplies by the appropriate factors 
of the AV theory: 





= A(A~ 


-1^ Hrpab 


f 


= A{A~ 






= A{A~ 





(4.2) 
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This form ensures that the transformation of both the current and its divergence satisfies 
the SUSY algebra. The consistency of this prescription with thermodynamics in the ideal 
case was a non-trivial test for its validity. 

Let us now comment on possible phenomenological implications of the results presented 
in this paper. The low energy theorems for theories with spontaneously broken symmetry 
tell us that the interaction between the goldstino and the supercurrent is suppressed by 
the SUSY breaking scale. However, in the paradigm of SUSY breaking in a hidden sector, 
it is well known that the effective interactions with the standard model particles should 
be suppressed by a higher scale, the mediation scale M, to ensure the decoupling of the 
goldstino in the M — > oo limit ||. In other words, when taking into account all quantum 
corrections, there should be some cancelations that cause the coupling of the gravitino 
to the SM plasma to be suppressed by -n. However, we suggest that this should be 
reconsidered. One may view the low energy fermionic zero mode as a linear combination of 
the hidden sector goldstino and the SM phonino. The hidden sector goldstino coupling is 
indeed suppressed by the mediation scale, but that of the phonino is not. Thus, although 
the contribution of the phonino to the wave function is smaller than that of the hidden 
sector goldstino, its interactions dominate and lead to the scattering terms discussed in 



section 3.2. Such enhanced interactions with the low energy modes of the plasma should 
be taken into account in the study of gravitino physics in the early universe, and perhaps 
even modify its production mechanism. To our knowledge, this effect was not considered 
in the literature to date. 

Another possible phenomenological consequence of this analysis might be in relaxing 
the constraints on the gravitino warm dark matter scenario: It is standard lore that dark 
matter cannot be a thermal gravitino relic - such a relic must have mass of the order of 
0.1 KeV to explain the dark matter mass density, while a thermal relic of the same mass 
would have a velocity distribution which is inconsistent with limits obtained from large 
scale structure [^lj. The lower bound on the mass of such a relic arising from large scale 
structure is of the order of several KeV, merely one order of magnitude higher than the 
mass suggested by the relic density. In fact, a borderline scenario of a KeV gravitino as a 
warm dark matter candidate has recently gained some renewed interest, since it provides a 
drop in the power spectrum which might be able to solve astrophysical inconsistencies such 
as the core vs. cusp question and the emptiness of voids (for a review we suggest p2|). 
One possible solution is to consider a 0.1 KeV gravitino which freezes out with the correct 
relic density, but with a velocity distribution which is modified by phonino interactions. In 
other words, when the gravitino is thermalized, one has to take into account the phonon- 
phonino effective interactions in the kinetic theory describing the freeze-out process of the 
gravitino. Since these interaction is proportional to the momentum of the phonon, this 
effect might slow down the phonino, and perhaps even lead to the required drop in the 
power spectrum. We leave this idea for future work. 
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A. Non-linearly realized supersymmetry 

In this appendix we review some well known results regarding the Akulov-Volkov (AV) 



formalism for non-linearly realized supersymmetry (NLRS) [12|, and give new general ex- 
pressions for conserved currents in this framework. In this section we will refer to the 
fermionic zero mode as a goldstino. We will assume a flat metric, use the conventions of 
1 14] and absorb the supersymmetry breaking scale into the definition of the goldstino. For 
additional reading we refer the reader to IS, 20]. 

A.l The standard non- linear realization of supersymmetry 

In the standard nonlinear realization of supersymmetry, the goldstino transforms as follows: 

5^G a = C + v^8^G a 

=f + v^G\ (A.l) 

where G is the goldstino and 

u£ = -i {Gcj^I - £a^G) . (A.2) 

Due to its non-linear transformation, the goldstino can be used to realize the SUSY trans- 
formation of other fields: 

5^ = v^d^ . (A.3) 

Indeed, using the transformation of 

6 v v» = -i (r,a^ - ^rj) + v»d u v» , (A.4) 

one can see that ( |A.3j ) satisfies the supersymmetry algebra 

(<$ € <$„ - <$„<5f) = -2i {ia^n ~ W^) . (A.5) 

In this discussion we will assume that there is no linear realization of supersymmetry. In 
other words, the energy scale is lower than the boson-fermion mass splitting, such that 
only one "representative" of a supermultiplet exists in the low-energy spectrum. 
A very useful object in this formalism is the AV vierbein: 

V = ^ + i { G ° ad »G ~ d^Ga a G) , (A.6) 

whose variation is 

^; = ^a;+4°H ( a - 7 ) 
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Its determinant A transforms as 



S^A = dJv^A\ , (Ai 



therefore the product of A with any objects which transforms according to the standard 
realization will transform into a total derivative. This is most useful when constructing 
actions invariant under NLRS. 



Gradients do not transform according to the standard realization ( A. 3 ). We therefore 
use the inverse vierbein, (A -1 )^ , which transforms as 

6s(A-\» = " (A-X»d v v» , (A.9) 

and define the AV covariant derivative 

V a cp = (A-V<V (A.10) 
which does transform according to the standard realization 

(V a 0) = V\3 V (V.tf) • (A.ll) 

The inverse vierbein itself can be written using the covariant derivative 

(A" 1 )/ = S^-ifGa^aG-VaGa^G) . (A.12) 
We note that the covariant derivatives are not commuting: 

[V M , V u ](f> = -i (v^Gcj t V v G - V„Ga T V M G) V t <P ■ (A.13) 
A. 2 Conservation equations in NLRS 

As explained in the text, conserved currents are unique objects in the sense that their 
supersymmetry transformation properties are more constrained - we require the trans- 
formation of both the current and its divergence to satisfy the SUSY algebra. Since the 
derivative in the conservation equations is not the AV covariant derivative, this is not a 
trivial condition. We therefore suggest the following form for the currents: 

f = A(A- 1 )/j a (A.14) 

where j a is the current defined such that it coincides with the non-SUSY current in the limit 
where the goldsitno is set to zero, and it transforms according to the standard realization. 
We refer to these currents as "standard realization" (SR) currents, defined in this 
way coincides with the classical current in the / — > oo limit as well, but it has different 
transformation properties: 

hf = d v {ylf -v»f) +v»d v f 
5 i (d^f) = (v»d u f) . (A.15) 
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Applying another NLRS transformation we find that the v dependent terms are symmetric 
under the exchange of £ and rj and will therefore vanish from the commutator: 



W = -* {v^Z ~ e°"v) duf + d p d v {vPvlf) - d v [vldpv^f) - d v (v»d p vlf 
We W 1 ) = -* - VI) (d u f) + d p d p (d v f)) ■ (A.i6) 

The important point is that with this construction the SUSY algebra is indeed satisfied, 
as in ( |A.5| ), for both the current and its divergence. 

In the next section of this appendix we will show that this is the form of the currents 
which is obtained from the Noether procedure. 

A. 3 Invariant actions and Noether currents in NLRS 

In this section we use the well-known construction of invariant actions in the AV formalism, 
compute Noether currents, and show that indeed they take the form suggested in ( |A.14| ). 
To make a Lagrangian invariant under non-linearly realized SUSY we make the following 
modification 

£(0 7 , dtf) -> ACty 1 , V<£ J ) (A.17) 

(the index / runs over all the fields in the theory) . In other words, we replace all derivatives 
with covariant derivatives and multiply the Lagrangian by A to make it transform into a 
total derivative as explained above . 

Let us first consider a Noether current corresponding to a global symmetry. We can de- 
fine a "standard realization" (SR) current applying the Noether formula to the Lagrangian: 



■a _ _V — Vj. )_L / A 18 

a(V a ^) da 

This is, of course, just a useful notation. The physical current is in fact given by 
3 ~ WM) da~- A[A )a diy^i) da— ~ A[A )aJ 

(A.19) 

in accordance with the prescription given in the previous section. 

We will now compute the stress-energy tensor. For this purpose we define an SR 
stress-energy tensor: 

* nh (A- 1 )^ d(AC) nb „ dC j , k . 

s - a d{{ A-v) = " c - mm vv (A - 20) 

In the limit of zero goldstino fluctuations it is identical to the canonical stress-energy 
tensor, and it transforms according to standard NLRS because it is constructed from 
Lorentz scalars and covariant derivatives: 

5^f ab = v^d p f ab . (A.21) 
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The canonical stress-energy tensor in the presence of goldstinos can be related to the 
classical one as follows: 

-In Ht A-Uap 9 (AC) 



-(A-X(A-±y 



A(A- l )ff ab (A.22) 



lb 

where we found the following relation useful: 

a ~ v c do^r^ • (A - 23) 

Again, we see that the form of the stress-energy tensor derived from the invariant action 
agrees with the prescription suggested above, therefore the transformation of both the ten- 
sor and its divergence satisfies the SUSY algebra. As in ( |A.15 ), the NLRS transformation 
is given by 

This tensor is not symmetric, but it can be made symmetric (on-shell) by adding a Belin- 
fante term 

Qiw — J'P^ _|_ ZQ [J^PP V _|_ flV"P _|_ JJVPP^ 

H pp, = € P^ T P aG(Tb G . (A.25) 

We move on to compute the supercurrent using the Noether procedure (as mentioned 
above, we assume that all fields transform according to the NLRS): 

& = gfeg) d(AC) a (ac) d&G ) a (ac) d(5^) _. Ar , „ n 

9£« 8(8^) d{d ^) d ^ W) d Z a lL - [aU) « 
= -2iT» a (a a G) a (A.26) 

Using the transformation properties of the stress-energy tensor and the goldstino we find 

W = % {v\S» - t£S£) + v»8 v S» - 2iT^ (a v l) a . (A.27) 

As required by the SUSY algebra, the supercurrent indeed transforms into the stress-energy 
tensor (plus derivative terms). Substituting ( A.22j ) into ( [A.26 ) we see that the supercurrent 
can also be written in terms of SR quantities dressed with appropriate goldstino factors. 

The conservation equations for the stress-energy tensor and the supercurrent combine 
to give the equation of motion for the goldstino: 

A(A~ 1 ) h ti f ab d^ (a a G) a = (A.28) 
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which reproduces the result appearing in ( 2.31 ). 



For completeness we also write the expression for the i?-current in cases in which only 
the goldstinos have R charge (e.g. hydrodynamics with NLRS, as discussed below): 

f R = -2T» a Ga a G 

k?R = d » {f^R ~ + v t d ^R ~ + ■ (A.29) 

The results of this section are consistent with [ 2C ] when using the goldstino self inter- 
actions Lagrangian AC = —f 2 A. 

B. Conservation of entropy current in fluids with phoninos 

As a consistency check for the identification of the entropy current in the supersymmetric 
case with the vector defined in ([^) , we show that the entropy is conserved by virtue 
of the thermodynamic relations and the equations of motion. We begin by writing the 
divergence of the entropy current as in ( |3.15| ) using the definition = A{A~ 1 )J J 'u a s and 
the thermodynamic relations ( |3.13| ) and fl3,14|) for the SR objects: 

-Td^ = -fy (A(A- X )*u a ) Ts - A{A- l ) a ^u a Td^s 

= -dp {A(A- r )f) (e + P)u a - A(A- r )f (e + P) d^u a - A(A~ 1 ) a fl u a d ti e 
= u a d, (A{A-\^ b ) - u a d, (A(A~ 1 ) a fl ) P . (B.l) 

Using the relation 

{A^A- 1 )*) = 2iA V bc (v a Ga c V b G - V h Ga c V a G) (B.2) 

we find 

—Tdfj, (W) = u a d^ (AiA-^f ^ - 2iu a AT bc (v a Ga c V b G - V b Ga c V a G^j (B.3) 

(Notice that the expression proportional to u a u b u c drops from this expression because the 
RHS in flB.2| ) is antisymmetric under exchange of a and b). Denoting p, = u a V a G 7 , and 
recalling ( [A.22j ) and ( [A. 26 ) we can rewrite this expression as: 



—Td^ {su") = u a d^ a - 2iT^ a (fia a d^G - d^Ga a fi) 

= (u a + 2% (fia a G - G<j a J) ) d^ a + fid^ + Jid^S 11 = (B.4) 

C. Phonon and phonino fluctuations 

In this appendix we collect some results we found useful for the computations of phonon 
and phonino interactions. 

The Lagrangian ( |3.28 ) by considering the following fluctuation in the matrix B IJ = 



sbu = b 2 /3 {dm j + dj^ + d^jd^j) (c.i) 



This notation is based on (2.27 
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The necessary derivatives of F are 



dF 



8B IJ 



lF B B l /3 8u 



8 2 F 



dB KL dB IJ 



j {F BB B 2 + F B B) (B- l ) LJ {B- l ) KL - ^F B B{B- x ) IK {B- x ) L j q 
\ {F BB Bl + F B B ) B- 4/3 8jj5 KL - ^F B B 1/3 5 IK 6 LJ (C.2) 



Upon integration by parts, a term djirjdjTTj becomes (f?/7Tj) 2 , giving the action ( [3.28 ). 

For the computation of the two-goldstinos one-phonon vertex, ( p. 33 ), we use the fol- 
lowing results: 



OBjj 



d[A°) 
dB 



d{A«) 
d 2 (AC) 



dA,?dB IJ 



-B 2 /3 [rja^j + r, aJ S») 
-B 5Z{l-5° a ) 

\f b B B q 2/3 (6Z6%J ~ c%j5Z (1 - S° a ) + (ri aI 6l + n aJ 5?)) 



(C.3) 



Computations of interaction terms with more phonons have an additional subtlety 
as they involve contributions from different orders in 5Bjj. For example, the goldstino 
interactions with two phonons has a contribution from the linear order in SBjj 



5 (AC) d -Ts ((Ga G - Ga G) + c 2 (Ga 1 ' djG - d I Ga I ~Gy) (d^jd^j) 
-iTs (Ga 1 djG - djG^G) d^md^nj , 

and another term from the 5Bjj5Bkl order 



(C.4) 



AC D -Ts (Jl + c 2 ) [GcjqG - Ga G) + d 2 (Ga 1 <9jG - d/GVG) J (djirj) 2 

+ l -Ts {[Gaot - Ga Q G) + c 2 s (Ga 1 ' d{G - d^a'G)) (djir K djTT K + djir K d K irj) 



-iTs (1 + c 2 ) (Ga'djG - djGa'G) (5/vrj + djin) (d K 7r K ) 
+ % -Ts (Ga'OjG - djGa'G) (d K iri + d^ K ) (d K 7Tj + djir K ) 

where we defined d 2 = (F BBB B 2 + 2F BB B) / F B and used 



(C.5) 



a 3 (ac) 



dA,?dB KL dB IJ 



F BB B 2 + F B B) 6ij6 kl - \f b B5 ik 5 L j) 



1 



B, 



-4/3 



€ (1 " <2 



FbbbB 6 + 2F BB B Z ) 6ijS kl - ]-F BB B 2 5 IK 5 L j 



+B 4/3 (F BB B 2 + F B B) Varf^KL + ^B va F B Bri aI 5^5 L j 



1 



-4/3 ■ 



(C.6) 



- 23 - 



References 

[1] A. K. Das and M. Kaku, Phys. Rev. D 18, 4540 (1978). 

[2] D. Boyanovsky, Phys. Rev. D 29 (1984) 743. 

[3] G. Policastro, JHEP 0902, 034 (2009) [arXiv:0812.0992 [hep-th]]. 

[4] J. P. Gauntlett, J. Sonncr and D. Waldram, JHEP 1111, 153 (2011) [arXiv: 1108. 1205 
[hep-th]]. 

[5] K. Kontoudi and G. Policastro, arXiv:1206.2067 [hep-th]. 
[6] R. G. Leigh and R. Rattazzi, Phys. Lett. B 352, 20 (1995) [hep-ph/9503402]. 
[7] K. Kratzcrt, Annals Phys. 308, 285 (2003) [hcp-th/0303260]. 
[8] B. dc Wit and D. Z. Freedman, Phys. Rev. Lett. 35, 827 (1975). 
B. do Wit and D. Z. Freedman, Phys. Rev. D 12, 2286 (1975). 

[9] K. Jensen, M. Kaminski, P. Kovtun, R. Meyer, A. Ritz and A. Yarom, arXiv:1203.3556 
[hep-th]. 

[10] N. Banerjee, J. Bhattacharya, S. Bhattacharyya, S. Jain, S. Minwalla and T. Sharma, 
arXiv:1203.3544 [hep-th]. 

[11] R. Loganayagam and P. Surowka, JHEP 1204, 097 (2012) [arXiv:1201.2812 [hep-th]]. 

[12] D. V. Volkov and V. P. Akulov, JETP Lett. 16, 438 (1972) [Pisma Zh. Eksp. Teor. Fiz. 16, 
621 (1972)]. 

[13] P. Kovtun and L. G. Yaffe, Phys. Rev. D 68, 025007 (2003) [hep-th/0303010]. 

[14] J. Wess and J. Bagger, Princeton, USA: Univ. Pr. (1992) 259 p 

[15] Z. Komargodski and N. Seiberg, JHEP 0909, 066 (2009) [arXiv:0907.2441 [hep-th]]. 

[16] D. T. Son, hep-ph/0204199. 

[17] S. Dubovsky, L. Hui, A. Nicolis and D. T. Son, arXiv:1107.0731 [hep-th]. 

[18] S. Endlich, A. Nicolis, R. Rattazzi, J. Wang, JHEP 1104, 102 (2011). [arXiv:1011.6396 
[hep-th]]. 

[19] T. E. Clark and S. T. Love, Phys. Rev. D 54, 5723 (1996) [hcp-ph/9608243]. 
[20] T. E. Clark and S. T. Love, Phys. Rev. D 39, 2391 (1989). 
[21] E. W. Kolb and M. S. Turner, Front. Phys. 69, 1 (1990). 

[22] H. J. de Vega, M. C. Falvella and N. G. Sanchez, arXiv: 1009.3494 [astro-ph.CO]. 



- 24 - 



